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Abstract. VIVA Hi observations of the Virgo spiral galaxy NGC 4501 are presented. The Hi disk is sharply 
truncated to the southwest, well within the stellar disk. A region of low surface-density gas, which is more 
extended than the main Hi disk, is discovered northeast of the galaxy center. These data are compared to existing 
6 cm polarized radio continuum emission. Ha, and optical broad band images. We observe a coincidence between 
the western Hi and polarized emission edges, on the one hand, and a faint Ha emission ridge, on the other. 
The polarized emission maxima are located within the gaps between the spiral arms and the faint Ha ridge. 
Based on the comparison of these observations with a sample of dynamical simulations with different values for 
maximum ram pressure and different inclination angles between the disk and the orbital plane, we conclude that 
ram pressure stripping can account for the main observed characteristics. NGC 4501 is stripped nearly edge-on, 
is heading southwest, and is ~ 200 — 300 Myr before peak ram pressure, i.e. its closest approach to M87. The 
southwestern ridge of enhanced gas surface density and enhanced polarized radio-continuum emission is due to 
ram pressure compression. It is argued that the faint western Ha emission ridge is induced by nearly edge-on ram 
pressure stripping. NGC 4501 represents an especially clear example of early stage ram pressure stripping of a 
large cluster-spiral galaxy. 

Key words. Galaxies: individual: NGC 4501 - Galaxies: interactions - Galaxies: ISM - Galaxies: kinematics and 
dynamics 



1. Introduction 

The Virgo cluster is dynamically young and spiral-rich 
making it an ideal laboratory for studying the influence 
of the cluster environment on spiral galaxies. Most of 
the Virgo spiral galaxies are Hi deficient, i.e. they have 
lost a significant amount of their ISM (Chamaraux et 
al. 1980, Giovanelli & Haynes 1983). Imaging Hi obser- 
vations have shown that these galaxies have truncated Hi 
disks (Cayatte et al. 1990). Thus, the cluster environment 
changes the Hi content and morphology of Virgo cluster 
spiral galaxies. However, it remains unclear at which clus- 
ter distance the spiral galaxies lose their external Hi disk. 
Based on deep Effelsberg Hi observations around 6 Virgo 
galaxies and a balance of previous detections of extrapla- 
nar gas in targeted galaxies VoUmer & Huchtmeier (2007) 
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propose a global picture where the outer gas disk (be- 
yond the optical radius R25) is removed much earlier than 
expected by the classical ram pressure criterion. They ar- 
gue that the vulnerable outer disk is removed much more 
rapidly than predicted by the Gunn & Gott criterion. 
Furthermore, Chung et al. (2007) found in a new Hi imag- 
ing survey of ^ 50 Virgo spiral galaxies 7 galaxies with 
long Hi tails pointing away from M87. These galaxies are 
located at intermediate distances from the cluster center 
(0.6-1 Mpc). They conclude that these one-sided Hi tail 
galaxies have recently arrived in the cluster, falling in on 
highly radial orbits. It appears that galaxies begin to lose 
their gas already at intermediate distances from the clus- 
ter center through ram-pressure or turbulent viscous strip- 
ping (Nulsen 1982) and tidal interactions with neighbours, 
or a combination of both. The outer gas disks of spiral 
galaxies at projected distances smaller than ~ 1 Mpc are 
thus affected by ram pressure. In this article we show that 
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this is also the case for the massive spiral galaxy NGC 4501 
(for its properties see table [1]). 

NGC 4501 is located at the lower end of the distance 
interval defined by Chung et al. (2007). Moreover, Cayatte 
et al. (1990, 1994) and Bohringer et al. (1997) proposed 
that it approaches the center of the Virgo cluster (M87) 
on a highly eccentric orbit. NGC 4501 is moderately Hi 
deficient and has an Hi disk whose southwestern side is 
truncated well within optical radius {R25 ) • The Hi surface 
density distribution within the disk is asymmetric with an 
overdense region to the southwest showing a sharp outer 
edge (Cayatte et al. 1990). Since the direction of these fea- 
tures is pointing toward the Virgo cluster center, Cayatte 
et al. (1990, 1994) and Bohringer et al. (1997) speculate 
that they are due to ram pressure compression (Gunn & 
Gott 1972). Since ram pressure only affects the interstellar 
medium of the galaxy and not its stars, one expects a sym- 
metric stellar disk. However, the DSS B band morphology 
is asymmetric, the southwestern side of the major axis 
being brighter than the northeastern side, rising the ques- 
tion if NGC 4501 is exclusively experiencing an interaction 
with the intracluster medium and if a past gravitational 
interaction can be excluded. 

In this article we investigate the nature of the inter- 
action of NGC 4501 with the cluster environment using 
VIVA (VLA Imaging of Virgo galaxies in Atomic gas, 
Chung et al. in prep.) Hi data, VLA 6 cm polarized ra- 
dio emission data, and a detailed dynamical model in- 
cluding MHD calculations. Whereas the Hi distribution 
and kinematics show a distorted gas disk where some of 
the gas is pushed to the north-east, the polarized radio- 
continuum emission represents an ideal tracer of ram pres- 
sure compression at the opposite side of the disk. The com- 
parison between observations and simulations shows that 
NGC 4501 is indeed approaching the cluster center. Ram 
pressure stripping can account for the main properties of 
this Virgo spiral galaxy. 

We describe the Hi and 6 cm polarized radio- 
continuum observations and their results in Sect. [5] and 
[3l Additional Ha data are discussed in Sect.[4l After a de- 
scription of the dynamical model, the best fitting model 
and following MHD calculations are presented (Sect.[5|). 
Our models are compared to the observational data in 
Sect. [5] followed by a discussion (Sect. [7]) and our conclu- 
sions (Sect. |8]). 

2. Observations 

2.1. VIVA Hi 21 cm line observations 

Hi line observations were made on January 20 1999, us- 
ing the C configuration of the VLA of the National Radio 
Astronomy Observatory (NRAOjll- The total on-source 
time was 5.6 hr. The bandpasses of 2x3.125 MHz were 
centered at 2280 km s~^ configured to produce two polar- 
izations and 63 channels with on-line Banning smoothing, 

^ NRAO is a facility of National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



Table 1. Physical Parameters of NGC 4501 



Other names M 88 

VCC 1401 
UGC 7675 

a (2000)^ 12''31™59.2= 

S (2000)^" 14°25'14" 

Morphological type" Sb 

Distance to the cluster center 2.0°, 0.6 Mpc 

Optical diameter Djg'' 6.9', 34 kpc 

B?," 9.86 

Systemic heliocentric velocity" (kms~^) 2280±4 

Velocity with respect to Virgo mean (kms-^) 1130 

Distance D (Mpc) 17 

Vrotmax (kms-^) 300*' 

PA 142° ^ 

Inclination angle 56° ^ 

HI deficiency" 0.56±0.2 



" RC3, de Vaucouleurs et al. (1991) 
^ Guharthakurta et al. (1988) 
° Cayatte et al. (1994) 



yielding a channel spacing of 10.4 km s^^. The nearby 
phase calibrator 1252+199 was observed every 25 minutes, 
and 1328-1-307 was used as a flux calibrator. 

We used the NRAO Astronomical Image Processing 
System (AIPS) for calibration and mapping. The contin- 
uum subtraction was done using a linear fit to the visibil- 
ity data of line- free channels. To make the image cube we 
applied a weighting scheme intermediate between uniform 
and natural, but closer to a natural weighting scheme (ro- 
bust=l), which resulted in a resolution of 17". To maxi- 
mize sensitivity, we made a second data cube with tapered 
data (resolution of 30"). The resulting rms are 0.7 mJy 
and 1.2 mJy per beam and per channel, respectively, which 
corresponds to 2.9/1.5 x 10^^ cm~^. 

2.2. 6 cm polarized radio-continuum observations 

NGC 4501 was observed for 3:55 h on December 3 2005 
with the Very Large Array (VLA) of the National Radio 
Astronomy Observatory (NRAO) in the D array configu- 
ration. The band passes were 2 x 50 MHz. We used 3C286 
as the flux calibrator and 1254-1-116 as the phase cali- 
brator, the latter of which was observed every 40 min. 
Maps were made for both wavelengths using the AIPS 
task IMAGR with R0BUST=3. The flnal cleaned maps 
were convolved to a beam size of 18" x 18". We ended up 
with an rms level of the linear polarization, taken to be the 
mean rms in Stokes Q and U, of 11 /iJy/beam. The distri- 
bution of 6 cm polarized continuum emission was already 
published in VoUmer et al. (2007). For each map point we 
determined the apparent polarization B-vector, defined as 
the E- vector rotated by 90° . Taking into account the pos- 
sible bias due to Faraday rotation (Wezgowiec et al. 2007) 
this represents a good approximation of the orientation of 
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the sky-projected mean regular magnetic field, accurate 
to about ±10°. 



3. Results 

The 30" resolution Hi distribution, velocity field, and ve- 
locity dispersion distribution are shown in Fig. [TJ The su- 
perposition of the Hi distribution on a B band and an 
H band (1.65 /ini) image are presented in Figs. [2] and [3l 
As observed by Cayatte et al. (1990) the Hi distribution 
shows a central hole and is truncated near the optical ra- 
dius i?25- The Hi surface density distribution within the 
disk is asymmetric with a high surface-density ridge in the 
southwest. The average Hi surface density in the north- 
eastern part of the disk is ~3 times lower than that of 
the southwestern ridge. We observe low surface-density 
Hi (^ 2 — 5 X 10^^ cm~^) to the northeast of the galaxy 
center. This gas was not detected by Cayatte et al. (1990). 

The velocity field of NGC 4501 is relatively symmetric 
parallel to the minor axis even in the low surface-density 
northeastern region. The rotation curve becomes flat at a 
smaller galactocentric radius in the southeast. The veloc- 
ity dispersion is symmetric with two maxima caused by 
beam smearing over the inner part of the galaxy with a 
steeply rising rotation curve. 

The comparison between the Hi surface density distri- 
bution and optical images (Fig. [2] and [3]) shows that the 
southwestern part of the Hi disk is sharply truncated to 
well within the stellar disk. The outermost Hi disk follows 
one isophotal contour in the southwest, but extends well 
beyond this contour in the northeast. This is also the case 
for a small (~ 20") straight low surface-density region in 
the north pointing to the northwest. As already mentioned 
in Sect, [l] NGC 4501 has an asymmetric surface bright- 
ness distribution in the B band with flocculent spiral arms 
showing pitch angles which are larger in the northeastern 
part of the galactic disk than in the southwestern part. In 
addition, the southwestern part of the disk is brighter than 
the northeastern part. However, the last B band isophot 
contour is symmetric. In addition, NGC 4501 has a sym- 
metric surface brightness distribution in the H band which 
shows the old stellar component. We thus conclude that 
the B band asymmetries are due to the star formation his- 
tory and do not reflect an asymmetry of the gravitational 
potential. 

The comparison between the 6 cm polarized radio- 
continuum emission distribution and the 17" resolution 
Hi surface density distribution is shown in Fig. 21 The 
polarized emission shows an enhanced ridge south of the 
galaxy center. This ridge coincides with the high column 
density Hi ridge and is located entirely within the Hi gas 
distribution. In addition, we observe a local maximum of 
polarized emission ~ 20" to the north of the galaxy cen- 
ter. In normal field spirals the polarized radio-continuum 
emission is found in interarm regions, with B-vectors par- 
allel to the adjacent spiral arms (Beck 2005). As this inner 
interarm region is not resolved in our radio polarization 
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Fig. 1. Hi observations of NGC 4501. Upper panel: surface 
density distribution. Contour levels are (1, 2, 3, 4, 5, 6, 
7, 8, 9, 10)x (0.125 Jy/beamkms-i or 1.5 x lO^o cm'^). 
The beam (30" x 30") is plotted in the lower left corner 
of the image. Middle panel: velocity field. Contours are 
from 2010 kms"^ (northwest) to 2590 kms~^ (southeast) 
in steps of 20 kms^^. Lower panel: velocity dispersion 
distribution. Contour levels are (10, 20, 30, 40, 50) kms^^. 
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Fig. 2. Greyscale: Hi surface density distribution (0 to 
1.25 Jy/beamkms~^; resolution: 30"). Contours: B band 
image from GOLD Mine (Gavazzi et al. 2003). 
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Fig. 3. Greyscale: Hi surface density distribution (0 to 
1.25 Jy/beamkms~^; resolution: 30"). Contours: H band 
image from GOLD Mine (Gavazzi et al. 2003). 

data, the origin of the polarized emission in the north re- 
mains unclear. The magnetic field vectors are compared 
to the optical B band images in Fig. O If the large-scale 
magnetic field is determined by the spiral arms, the mag- 
netic field vectors should follow the spiral structure. This 
is what we observe in the inner ~ 1' = 5 kpc. However, 
the magnetic field vectors cross the spiral arm in the most 
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Fig. 4. Greyscale: Hi surface density distribution (0 to 
0.65 Jy/beamkms^^; resolution: 17"). Contours: 6 cm po- 
larized radio-continuum emission (VoUmer et al. 2007). 
The contour levels are (0.7, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
10) X 72 /iJy/beam. The vectors of the magnetic field 
are uncorrected for Faraday rotation. Their size is pro- 
portional to the intensity of the polarized emission. The 
beam (18" x 18") is plotted in the lower left corner of the 
image. 

northeastern part of the polarized emission distribution. 
This is similar to the behavior of the large-scale mag- 
netic field in NGC 3627 (Soida et al. 2001) where a mag- 
netic arm crosses an optical arm, which is interpreted as 
a sign of tidal interactions within the Leo Triplet group. 
On the other hand, the magnetic field within the ridge of 
enhanced polarized radio-continuum emission follows the 
spiral structure. Fig. [S] shows profiles parallel to the mi- 
nor axis centered on the major axis 1' southeast of the 
galaxy center (see Fig. [5]). We show the Hi surface density, 
the 6 cm polarized radio-continuum emission, and the Ha 
emission. The Hi and the polarized emission profiles rise 
together and peak at the same location. The polarized 
radio-continuum emission is enhanced by ram pressure 
compression of the magnetic field. If the original magnetic 
field was random, the compressed field is anisotropic and 
we do not expect Faraday rotation. A detailed analysis 
of Faraday rotation measures of Virgo Cluster spirals will 
be subject of future studies. The intensity ratio between 
the peak and the base of the polarized emission can be 
used to measure the compression ratio (Beck et al. 2005). 
The high resolution faint Ha emission coincides with the 
Hi and polarized emission peaks. A similar behavior is 
also observed in the interacting spiral galaxy NGC 2276 
(Hummel & Beck 1995) where a ridge of enhanced polar- 
ized radio-continuum emission is observed in the south of 
the galaxy, in the group spiral galaxy NGC 2442 (Harnett 
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Fig. 5. Greyscale: B band image from GOLD Mine 
(Gavazzi et al. 2003). The vectors of the magnetic field are 
uncorrected for Faraday rotation (resolution: 18"). Their 
size is proportional to the intensity of the polarized emis- 
sion as in Fig. |4l The line indicates the slit position for 
the profile of Fig. [6l 
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Fig. 6. Profiles parallel to the minor axis ^ 1' southeast 
of the galaxy center (see Fig. [5|). Solid: Hi surface density 
(resolution: 30"); dash-double-dotted: Hi surface density 
(resolution: 17"); dotted: 6 cm polarized radio-continuum 
emission (resolution: 18"); dashed: Ha emission at a res- 
olution of 30"; dash-dotted: Ha emission at a 2" resolu- 
tion. The y-axis has arbitrary units. Left corresponds to 
the southwest. The position of the major axis is marked 
with the vertical line. 



et al. 2004) where the northern polarized emission ridge is 
shifted towards the outer edge of the star-forming region, 
and the Virgo spiral galaxy NGC 4254 (Soida et al. 1996) 
where the Ha-bright spiral arm is placed inwards with 



Fig. 7. Greyscale: Hi surface density distribution (0 to 
0.65 Jy/beam kms^^; 17" resolution). Contours: Ha emis- 
sion distribution from GOLD Mine (Gavazzi et al. 2003). 
The contour levels are chosen in a way to make the faint 
southwestern Ha arm visible. 



respect to the compressed region. In NGC 2276 a tidal 
compression due to a galaxy-galaxy interaction is respon- 
sible for the enhanced polarized emission. For NGC 2442 
and NGC 4254 a tidal interaction also cannot be excluded. 



4. Star formation 

We use an Ha image to compare the Hi and polarized 
radio-continuum emission to the distribution of massive 
star formation (Figs. [7] and [8|). The Ha distribution shows 
a ring structure in the inner disk (R < 40") and a promi- 
nent spiral arm in the northwest. The spiral arm begins 
north of the galaxy center, continues to the northwest, 
bends to the southwest, and finally follows the southwest- 
ern edge of the Hi distribution. There is the beginning of a 
second Ha spiral arm in the southeast of the galactic disk. 
This second spiral arm abruptly ends ~ 1' = 5 kpc from 
the galaxy center. The comparison between the Ha emis- 
sion distribution and the 6 cm polarized radio-continuum 
emission distribution (Fig. [8]) shows that the maximum of 
the polarized emission south of the galaxy center is located 
between the inner Ha ring and the outer Ha arm. The lo- 
cation of maximum polarized radio-continuum emission 
may indicate the region of strongest field compression, co- 
inciding with maximum HI emission. Alternatively, the 
field could be tangled in star forming regions located in 
the Ha arms. Note that the limited resolution of the radio 
observations leads to beam depolarization in the inner disk 
around the major axis if the field lines follow the strong 
curvature of the optical spiral structure. The second lo- 
cal maximum northwest from the first maximum is also 
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Fig. 8. Greyscale: Ha emission distribution. Contours: 
6 cm polarized radio-continuum emission. The contour 
levels are (1, 2, 4, 6, 8, 10) x 72 /xJy/beam. 

located in an Ha emission gap. On the other hand, the 
two local maxima located at ~ 20" north and south of the 
galaxy center do not coincide with such clear gaps. The 
maxima in the ridge of enhanced polarized emission are 
thus different from those of the inner disk, which are pos- 
sibly due to the shear caused by the spiral arms. Spitzer 
8 /Ltm and 24 /xm maps (in prep.) show essentially the same 
spiral arm features as the Ha. 

5. Models 

To interpret the observations described in Sect.[3l we com- 
pare them to a detailed dynamical model which includes 
ram pressure. In a second step, we solve the induction 
equation for the velocity fields of the dynamical model 
to determine the evolution of the large-scale magnetic 
field. Assuming a distribution of relativistic electrons, we 
then calculate the model distribution of polarized radio- 
continuum emission. 

5.1. Dynamical model 

We use the N-body code described in Vollmer et al. (2001) 
which consists of two components: a non-coUisional com- 
ponent that simulates the stellar bulge/disk and the dark 
halo, and a coUisional component that simulates the ISM. 

The non-coUisional component consists of 49 125 par- 
ticles, which simulate the galactic halo, bulge, and disk. 
The characteristics of the different galactic components 
are shown in Table [21 The resulting rotation velocity is 
~300 kms-^ 

We have adopted a model where the ISM is simulated 
as a coUisional component, i.e. as discrete particles which 



Table 2. Total mass, number of particles N, particle mass 
M, and smoothing length I for the different galactic com- 
ponents. 



component 


Mtot (Mq) 


N M (Mq) / (pc) 


halo 


7.310" 


16384 4.5 10' 1200 


bulge 


2.510^° 


16384 1.5 10** 180 


disk 


1.310" 


16384 7.910^ 240 



possess a mass and a radius and which can have inelas- 
tic collisions (sticky particles). Since the ISM is a turbu- 
lent and fractal medium (see e.g. Elmegreen & Falgarone 
1996), it is neither continuous nor discrete. The volume 
filling factor of the warm and cold phases is smaller than 
one. The warm neutral and ionized gas fill about 30 — 50% 
of the volume, whereas cold neutral gas has a volume fill- 
ing factor smaller than 10% (Boulares & Cox 1990). It 
is not clear how this fraction changes, when an external 
pressure is applied. In contrast to smoothed particles hy- 
drodynamics (SPH), which is a quasi continuous approach 
and where the particles cannot penetrate each other, our 
approach allows a finite penetration length, which is given 
by the mass-radius relation of the particles. Both methods 
have their advantages and their limits. The advantage of 
our approach is that ram pressure can be included eas- 
ily as an additional acceleration on particles that are not 
protected by other particles (see Vollmer et al. 2001). In 
this way we avoid the problem of treating the huge density 
contrast between the ICM (n ^ 10^^ cm"^) and the ISM 
(n > 1 cm~'^) of the galaxy. 

The 20 000 particles of the coUisional component repre- 
sent gas cloud complexes which are evolving in the grav- 
itational potential of the galaxy. The total assumed gas 
mass is M*"* = 8.9 10^ M©, which corresponds to the to- 
tal neutral gas mass before stripping. Following Vollmer & 
Huchtmeier (2007) we assume an initial gas disk which is 
truncated at approximately the optical radius (~ 17 kpc). 
To each particle a radius is attributed depending on its 
mass. During the disk evolution the particles can have in- 
elastic collisions, the outcome of which (coalescence, mass 
exchange, or fragmentation) is simplified following Wiegel 
(1994). This results in an effective gas viscosity in the disk. 

As the galaxy moves through the ICM, its gas clouds 
are accelerated by ram pressure. Within the galaxy's iner- 
tial system its clouds are exposed to a wind coming from a 
direction opposite to that of the galaxy's motion through 
a static ICM. The temporal ram pressure profile has the 
form of a Lorentzian, which is a good approximation for 
galaxies on highly eccentric orbits within the Virgo cluster 
(Vollmer et al. 2001). The effect of ram pressure on the 
clouds is simulated by an additional force on the clouds in 
the wind direction. Only clouds which are not protected 
by other clouds against the wind are affected. 

The particle trajectories are integrated using an adap- 
tive timestep for each particle. This method is described 



VoUmer et al.: Pre-peak ram pressure stripping in NGC 4501 



7 




time (Gyr) 



galactic plane 




wind direction 



Fig. 9. Galaxy geometry. Upper panel: image plane. PA: 
position angle; incl: inclination of the galactic disk with 
respect to the plane of the sky; a: azimuthal viewing angle 
within the galaxy plane defined by PA and incl. Lower 
panel: perpendicular view from the top. i: disk-wind angle. 

in Springel et al. (2001). The following criterion for an 
individual timestep is applied: 



20 kms" 



(1) 



where is the acceleration of the particle i. The mini- 
mum value of ti defines the global timestep used for the 
Burlisch-Stoer integrator that integrates the coUisional 
component. 

5.2. Search for the best fit model 

In this section we will constrain the parameters of the ram 
pressure stripping event which are (i) the peak ram pres- 
sure, (ii) the temporal ram pressure profile, (iii) time since 
peak ram pressure, (iv) the disk-wind angle between the 
galaxy's disk and the intracluster medium wind direction, 
and (v) the azimuthal viewing angle for the observed in- 
clination and position angles (Fig. [5]). These parameters 
are related to the observed quantities which are (1) the 



Fig. 10. The adopted temporal ram pressure profiles for 
the low ram pressure (LRP) and the high ram pressure 
(HRP) model. Thick lines: ram pressure profiles with dou- 
bled widths (LRPl/2, HRPl/2). The vertical hues show 
the adopted timesteps for the comparison between model 
and observations. 

position angle PA, (2) the inclination angle of the galac- 
tic disk incl, (3) the line-of-sight velocity of the galaxy 
with respect to the cluster mean, and (4) the projected 
ICM wind direction. The position angle and inclination of 
NGC 4501 define the galaxy plane in three dimensional 
space. The model galaxy can then be rotated within this 
plane by the azimuthal viewing angle (see below). The 
three dimensional model wind direction, the line-of-sight 
velocity of the galaxy, and the projected ICM wind direc- 
tion are thus functions of the azimuthal viewing angle. 

We assume a highly eccentric galaxy orbit within the 
Virgo cluster leading to a temporal ram pressure profile 
of the form: 



Pr 



Pn 



■-HW 



t^ + fi 



(2) 



HW 



where tnw is the width of the profile (VoUmer et al. 2001). 
We define t = Myr as the time when ram pressure is 
maximum. As shown in Vollmer et al. (2001) the width of 
the profile depends on the ram pressure maximum. For a 
high maximum ram pressure the galaxy has to come very 
close to the cluster center and consequently its velocity 
is higher while it crosses the high IGM density region in 
the cluster core. We use two typical profiles for strong 
and moderately strong ram pressure stripping (Fig. llOi 
see Vollmer et al. 2001): 



80 Myr and 



1. low ram pressure model (LRP): 
Pmax = 2000 cm-3(kms"l)^ tnw 

2. high ram pressure model (HRP): 

= 5000 cm-3(kms-^)^ t^w = 50 Myr. 

It turned out that, in both simulations, NGC 4501 en- 
counters peak ram pressure within 120 Myr. Since this is 
not consistent with its location within the cluster assum- 
ing classical galaxy orbits in a static smooth intracluster 
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Fig. 11. Left panel: galaxy orbit in the Virgo cluster 
leading to temporal ram pressure profiles with a larger 
width. Right panel: resulting ram pressure profiles (solid 
lines) and LRPl/2, HRPl/2 ram pressure profiles (dashed 
lines) . 

medium (Vollmer et al. 2001), we made a second set of 
simulations for which we doubled the widths of the tem- 
poral ram pressure profiles fFig. llOp : 

1. low ram pressure model (LRPl/2): 

p^^^ = 2000 cm-3(kms"i)2, tnw = 160 Myr and 

2. high ram pressure model (HRPl/2): 

p^^^ = 5000 cm-3(kms-i)2, tnw = 100 Myr. 

In these simulations the galaxy does not go as close to 
the cluster core (Fig. [TT|) as in the LRP and HRP mod- 
els. The minimum distance to the cluster center is about 
10% larger than that of the orbits shown in Vollmer et 
al. (2001). The minimum distance of the LRPl/2 is about 
0.25 Mpc, that of the HRPl/2 model is about 0.1 Mpc. 
The larger width of the ram pressure profile is mainly due 
to a broader intracluster medium distribution compared 
to the LRP and HRP simulations. 

With these 4 simulations we aim to reproduce the main 
characteristics of our Hi observations: 

— the low surface-density gas located northeast of the 
galaxy center, 

— the region of high surface-density gas at the opposite 
windward side, 

— the asymmetry along the major axis with a more ex- 
tended emission to the southeast, 

— the rising rotation curve to the northwest, 

— the southwestern ridge of polarized radio-continuum 
emission. 

The first snapshot selection is mainly based on the Hi 
morphology. For the LRP and HRP simulations we iden- 
tify one best fit model (LRP: Fig. [13l HRP: Fig.[T4l), for 
the LRPl/2 and HRPl/2 simulations two best fit models 
(LRPl and LRP2: Fig.[l5l HRPl and HRP2: Fig. [HI) are 
selected. In a second step (Sect. [6]), gas kinematics and the 
polarized radio contiuum emission distribution are taken 
into account. Table [3] summarizes our findings. 

Ram pressure efficiency depends on the disk-wind an- 
gle i between the galactic disk and the ICM wind direction 
(Vollmer et al. 2001). A disk-wind angle of z < —30° is not 
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Fig. 12. Components of the 3D direction of the galaxy's 
velocity as a function of the azimuthal viewing angle a. 
The length of the vector is unity. Solid line: right as- 
cension; dotted line: declination; dashed line: line-of-sight 
direction for four different disk-wind angle between the 
galactic disk and the ICM wind i = -30°, i = -10°, 
i = 10°, and i = 30°. The vertical line represents a view- 
ing angle of 140°. 

compatible with NGC 450 1's positive radial velocity with 
respect to the Virgo cluster. Disk-wind angles i > 30° lead 
to projected wind directions to the northeast for almost 
all azimuthal viewing angles. However, the Hi and po- 
larized emission observations indicate that the projected 
wind direction has an important component to the west. 
Therefore, we do not consider stripping with |z| > 30°. For 
each model (LRP and HRP) we made 4 simulations with 
4 different disk-wind angles between the galaxy's disk and 
the ICM wind direction: (i) i = -30°, (ii) i = -10°, (iii) 
i — 10°, and (iv) i ~ 30°. A disk- wind of i = 0° means 
that the galactic disk is parallel to the ICM wind direc- 
tion. The azimuthal viewing angle, which is the angle in 
the galactic plane defined by the galaxy's inclination and 
position angles, is chosen in a way to fit the observed Hi 
distribution, i.e. the position of the extended northeast- 
ern low surface-density gas and to reproduce the positive 
line-of-sight component of the wind direction (the galaxy 
is moving away from the observer). In Figs. [12] we show 
the components of the 3D ICM wind direction. For the 
chosen azimuthal viewing angle a = 140° the line-of-sight 
velocity is positive and represents between 50% and 100% 
of the total velocity with respect to the cluster mean. 

In Figs . [T51 and [Til we show snapshots of the LRP and 
HRP model at 3 different timesteps -150 Myr, -100 Myr, 
and —50 Myr, i.e. before the ram pressure maximum. 
All model gas distributions a,t t = —150 Myr have sym- 
metric outer contours and a local surface-density maxi- 
mum in the southeast, which is due to the internal spi- 
ral structure. At this timestep ram pressure is not strong 
enough to affect the gas distribution of the disk signif- 
icantly. However, gas of lower surface density is already 
displaced for i < -150 Myr. At i = -120 Myr the south- 



VoUmer et al.: Pre-peak ram pressure stripping in NGC 4501 
Table 3. Features which are reproduced by the different models. 



9 



model 


LRP 


HRP 


LRPl 


LRP2 


HRPl 


HRP2 


NE low surface density gas 


+ 


+ 


+ 


+ 


+ 


+ 


SW high surface density gas ridge 


+ 


+ 


+ 


+ 


+ 


+ 


HI asymmetry along major axis 








+ 




+ 


NW rising rotation curve 


+ 


+ 


+ 


+ 


+ 


+ 


SW ridge of polarized radio emission 


+ 


+ 


+ 


+ 


+ 


+ 


symmetric velocity dispersion distribution 


+ 


+ 


+ 




+ 




consistency with galaxy orbit 






+ 


? 


+ 


? 



western edge of the gas disk is now compressed and shows 
a sharp edge. On the opposite side, i.e. in the northeast, 
gas is pushed to larger galactic radii. Since the gas is still 
rotating there, it forms a kind of detached low surface spi- 
ral arm. This arm is located somewhat closer to the disk 
for i — 30° than for i — —30°, because it is located above 
the disk which makes its projected distance smaller. At 
t — —90 Myr the northeastern low surface arm has moved 
to the southeast and its projected distance from the main 
disk has increased. The LRP and HRP models are indistin- 
guishable at t = -150 Myr. At i = -120 Myr the overall 
Hi morphologies of the LRP and HRP model snapshots 
are similar and it is mainly the structure of the surface 
density distribution of the northeastern arm and the gas 
disk which change from one model to another. The latter 
statement also applies for the snapshots at < = —90 Myr. 
Based on these snapshots, only a simulation snapshot at 
t = —120 Myr reproduces the observed southwestern com- 
pression region and the northeastern low surface-density 
arm. In all model snapshots at i = — 120 Myr we observe 
a region of high surface density ^ 5 kpc southwest of the 
galaxy center, which is not present in the Hi data. Since 
this is due to the internal structure of the model galaxy, 
which (i) is not affected by ram pressure because of its 
location in the inner galactic disk and (ii) depends on the 
initial conditions, we do not take it into account for the 
comparison between models and observations. Based on 
a detailed comparison of the model gas distribution, in- 
cluding the southwestern overdensities and the northeast- 
ern low surface-density arm, we chose two final snapshots 
with i = 10° for the LRP model and i = -10° for the 
HRP model labeled in Figs. [13] and [H 

Both simulations yield a timescale of 120 Myr be- 
fore peak stripping. The galaxy is moving to the south- 
west with a velocity component in the sky-plane of 
1000 kms~^. Until peak ram pressure it will cover a 
distance of 0.13 Mpc. Its actual distance from the clus- 
ter center at peak ram pressure will be greater than 
0.5 Mpc. However, typical orbits associated to the ram 
pressure profiles LRP and HRP have impact parameters of 
0.1-0.2 Mpc (VoUmer et al. 2001). To reconcile this differ- 
ence, we doubled the widths of the temporal ram pressure 
profiles (LRP 1/2 and HRP 1/2 simulations; Fig.IHl)- This 
is still within the uncertainty of the intracluster medium 



density distribution given by Schindler et al. (1999). The 
impact parameters for orbits associated with these ram 
pressure profiles are 0.2 — 0.4 Mpc. Since i = 30° lead to 
negligible velocity components in the plane of the sky for 
the adopted azimuthal viewing angle fFig. [T2|) . we do not 
consider this possibility anymore. 

In Figs. [TSl and [TBI we show snapshots of the LRPl/2 
and HRPl/2 model at 5 different timesteps —380 Myr, 
-320 Myr, -260 Myr, -190 Myr, and -130 Myr, i.e. be- 
fore ram pressure maximum. The evolution of the LRP 1/2 
and HRP 1/2 models for a given disk- wind angle between 
the disk and the ram pressure wind are qualitatively 
similar. For timesteps t < —150 Myr the Hi distribu- 
tion of the different disk-wind angles are also similar. At 
t = —380 Myr there is more high surface density in the 
northwestern part of the galaxy, whereas we observe a low 
surface-density arm in the southeast. At t = —320 Myr 
the disk emission shows is approximately symmetric along 
the major axis and a low surface-density arm is formed 
east of the galaxy center. This snapshot corresponds to 
the LRP and HRP snapshots at t ^ -120 Myr (Figs. [H 
and I14p . because the ram pressure induced momentum 
transfer at the two respective timesteps is about the same 
(Fig. [T0|) . At t = —260 Myr this arm is located closer to 
the galaxy center and again a southeastern low surface- 
density arm is formed. A pronounced asymmetry along 
the minor axis is observed at t — —190 Myr with more 
emission in the compressed western part of the disk. In 
addition, we find a less pronounced asymmetry along the 
major axis with more emission in the southeastern part. 
At t = —130 Myr the gas in the western part of the disk 
is pushed to smaller radii and we observe a low surface- 
density arm in the southeast of the galaxy center. Based on 
this set of simulations, 4 snapshots with i = —10° repro- 
duce our Hi observations: LRPl - t = -320 Myr, LRP2 
- t ^ -190 Myr, HRPl - t = -320 Myr, and HRP2 - 
t — —190 Myr. All 4 snapshots show a low surface-density 
gas located northeast of the galaxy center and an region 
of high surface-density gas at the opposite windward side. 
Only the snapshots at t — —190 Myr show an asymmetry 
along the major axis with a more extended emission to 
the southeast. 

The time evolutions of these simulations are shown 
in Fig. 1171 Since ram pressure only affects the interstel- 
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Fig. 13. Model snapshots of the gas column density for different simulations of the low ram pressure (LRP) model. 
Contours are (1, 2, 3, 4, 5, 6, 7, 8, 9, 10) x4 10^" cm^^. The ram pressure profile for all simulations isq given by Eq. [2l 
The disk-wind angle i between the orbital and the galaxy's disk plane and the timestep of the snapshot are varied. 
Left column: t = —150 Myr. Middle column: t = —100 Myr. Right column: t — —50 Myr. Upper row: i ~ —30°. 
Middle rows: i = —10° and i = 10°. Lower row: i = 30°. The arrow represents the projected ICM wind direction. The 
absence of the arrow for i — 30° is due to a zero projected galaxy velocity. 



lar medium of the galaxy, the stellar disk stays symmet- 
ric during the whole simulation. The sky-projected wind 
direction is northeastward for the LRPl/2 and HRPl/2 
model. As expected, the stripping radius of the HRPl/2 
model is smaller than that of the LRPl/2 model, i.e. the 
galaxy in the HRPl/2 model loses more gas than that of 
the LRPl/2 model. 



5.3. Magnetic field evolution 

Otmianowska-Mazur & VoUmer (2003) studied the evolu- 
tion of the large-scale magnetic field during a ram pres- 
sure stripping event. They calculated the magnetic field 
structure by solving the induction equation on the veloc- 
ity fields produced by the dynamical model. The polarized 
radio-continuum emission has been calculated by assum- 
ing a Gaussian spatial distribution of relativistic electrons. 
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Fig. 14. Model snapshots of the gas column density for different simulations of the high ram pressure (HRP) model. 
Contours are the same as in Fig. [131 The ram pressure profile for all simulations is given by Eq. [2] The disk-wind 
angle i between the orbital and the galaxy's disk plane and the timestep of the snapshot are varied. Left column: 
t = -150 Myr. Middle column: t = -100 Myr. Right column: t = -50 Myr. Upper row: i = -30°. Middle rows: 
i = —10° and i — 10°. Lower row: i = 30°. The absence of the arrow for i — 30° is due to a zero projected galaxy 
velocity. 



This procedure allowed them to study the evolution of the 
observable polarized radio continuum emission during a 
ram pressure stripping event. 

We apply the same procedure as Otmianowska-Mazur 
& VoUmer (2003) to a similar ram pressure stripping event 
(Sect. [5|). The Zeus3D code (Stone & Norman 1992a and 
b) is used to solve the induction equation: 

dB/dt ^Tot{v X B) -Toti-q TotB) (3) 



where B is the magnetic induction, v is the large-scale ve- 
locity of the gas, and 77 is the coefficient of a turbulent dif- 
fusion. We use a physical diffusion of ry = 5 x 10^^ cm^s"^ 
(Elstner et al. 2000). We do not implement any dynamo 
process. The initial magnetic field is purely toroidal with 
a strength of 10 fiG. 

The induction equation is solved on rectangular co- 
ordinates (XYZ). The number of grid points used is 
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Fig. 15. Model snapshots of the gas column density for different simulations of the low ram pressure (LRPl/2) model 
with the large temporal ram pressure profile. Contours are the same as in Fig. 1131 The ram pressure profile for all 
simulations is given by Eq. [51 The disk-wind angle i between the orbital and the galaxy's disk plane and the timestep 
of the snapshot are varied. 
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Fig. 16. Model snapshots of the gas column density for different simulations of the high ram pressure (HRPl/2) model 
with the large temporal ram pressure profile. Contours are the same as in Fig. 1131 The ram pressure profile for all 
simulations is given by Eq. [51 The disk-wind angle i between the orbital and the galaxy's disk plane and the timestep 
of the snapshot are varied. 
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Fig. 17. Evolution of the model stellar (1st column) and 
gas disk (2nd column). The contours of the stellar distri- 
bution are logarithmic, those of the gas distribution are 
the same as in Fig. [T31 The arrow indicates the direction 
of ram pressure, i.e. it is opposite to the galaxy's veloc- 
ity vector, and its size is proportional to pWgai- Maximum 
ram pressure occurs at i = Myr. The timestep of each 
snapshot is marked in each panel showing the stellar disk. 



215x215x91 along the X, Y and Z axis, respectively with 
a grid size of 400 pc, resulting in a size of the modeled box 
of 86 kpc X 86 kpc x 36 kpc. Since the N-body code is 
discrete whereas the MHD code is using a grid, we have 
to interpolate the discrete velocities on the grid. This is 
done using a method known as "Kriging" with a density- 
dependent smoothing length (Isaaks & Srivastava 1989). 
It turned out that we had to use a large smoothing length 
to suppress the noise in the velocity field of the outer disk, 
which is due to a small local particle density. In this way 
we avoid numerical artifacts in the velocity field at the 
outer disk (see Soida et al. 2006) . As the spline interpola- 
tion used in Otmianowska-Mazur & VoUmer (2003), this 
has the consequence that the velocity field at the edge of 
the gas distribution is more extended than the gas distri- 
bution itself. Since there are gradients in this velocity field 
due to differential rotation, induction leads to a magnetic 
field which extends beyond the edge of the gas distribu- 



tion. This affects the polarized emission beyond the gas 
distribution, but not inside, which is what we are inter- 
ested in. 

The evolution of the polarized radio-continuum emis- 
sion without Faraday rotation is presented in Fig. 1181 
The timesteps are the same as in Fig. [TT] The disk ro- 
tates counter-clockwise. The gas surface density, which 
is smoothed to a resolution of ^ 100 pc, is shown in 
greyscales, the polarized radio-continuum emission as con- 
tours, and the magnetic field vectors projected on the 
plane of the sky as lines. We assume a Gaussian distribu- 
tion of relativistic electrons in R and z directions: nrci = 
riQ exp ( — (r/rpi)^) exp ( — (z/r^)^), where = 6.3 kpc 
and Tz = 1.0 kpc. This translates into a FWHM of 5.2 kpc 
and 0.83 kpc, respectively. By assuming this smooth distri- 
bution we imply no local equipartition between the energy 
densities of total cosmic rays and total magnetic field. The 
galaxy moves to the lower right corner, i.e. the ram pres- 
sure wind comes from this direction. Our grid size does not 
permit us to resolve the rotation in the inner part of the 
galaxy. The errors due to discretization lead to an artifi- 
cial radial diffusion of the magnetic field out of the galactic 
disk. During the evolution without a significant external 
influence (-380 Myr< t < -290 Myr) the model distri- 
bution of polarized radio-continuum emission is approxi- 
mately symmetric. Once ram pressure compresses the gas, 
an enhanced ridge of polarized emission appears in the 
southwest. At t = —230 Myr the internal magnetic field is 
lost due to the artificial diffusion, and only the compres- 
sion induced ridge of polarized emission is visible. The 
evolution of polarized radio-continuum emission is very 
similar in the LRPl/2 and HRPl/2 models. 



6. Comparison between the models and 
observations 

In this section we compare the observed properties of 
the gas, star formation, and large-scale magnetic field of 
NGC 4501 with the properties of the simulated galax- 
ies. Assuming a static intracluster medium, the projected 
direction of the ram pressure wind indicates the pro- 
jected direction of the galaxy's motion within the intr- 
acluster medium. For the relevant simulations IIRPl/2 
and LRPl/2 with i = —10° this direction is southwest 
fFigs.[T51and [T5)l . A time to peak ram pressure of 320 Myr 
or 190 Myr leads the galaxy to a distance of 0.4 Mpc 
and 0.5 Mpc from the cluster center at peak ram pres- 
sure assuming a linear orbit with a constant velocity (see 
Sect. [7]). Projection effects will increase, orbital curva- 
ture and acceleration will decrease these values. An in- 
spection of simulated galaxy orbits (Fig. [TT|) shows that 
the difference between the location at peak ram pressure 
and the extrapolated location assuming a linear orbit is 
0.1 Mpc for a timescale of 320 Myr (~ 0.03 Mpc for 
a timescale of 190 Myr). Our adopted minimum distance 
scale is thus 0.3-0.4 Mpc. This distance is only compat- 
ible with the low ram pressure simulations LRPl/2 (see 
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Fig. 18. Evolution of the polarized radio-continuum emission. Time steps are indicated in the upper left corner of 
each panel. The disk is seen with the position and inclination of NGC 4501 (Table [T]) and rotates counter-clockwise. 
Greyscale: gas surface density in logarithmic scale. Contours: polarized radio-continuum emission. The magnetic field 
vectors are superimposed onto the gas surface density. The assumed beamsize for the polarized radio-continuum 
emission is 20". 



Sect. 15. 2( ). Therefore, in the following we will only compare 
the LRPl/2 model to our multi- wavelength observations. 

6.1. The gas moment maps 

As described in Sect. [3] the observed Hi distribution fol- 
lows one isophot contour of the optical image except in 
the northeast and the north. Our model snapshots LRPl 
and LRP2 show the same qualitative behavior (Fig. [T9|) . 
However, the LRP2 model gas distribution in the west is 
less extended than the observed Hi distribution, i.e. the 



model gas is pushed to smaller galactic radii than it is ob- 
served. On the other hand, the observed asymmetry along 
the major axis is only reproduced by the LRP2 model. 

Fig. [20] compares the gas distribution, velocity field, 
and the distribution of the velocity dispersion of our Hi 
observations with our simulations. Both models repro- 
duce qualitatively the overall observed Hi surface den- 
sity distribution of the outer gas disk. Both models show 
an overdensity in the compressed southwestern region of 
the disk. Both models also show an extended low surface- 
density region in northeast, as it is observed. There are, 



16 



VoUmer et al.: Pre-peak ram pressure stripping in NGC 4501 






1biiiFbm"iT1^ 



Fig. 19. Comparison between the gas surface density distribution (greyscale) and the stellar surface density distribution 
(contours). Left panel: B band contours on Hi surface density distribution (resolution: 30"). In all figures the greyscales 
use the full dynamic ranges of the gas distribution. 



however, quantitative differences: the northeastern low 
surface-density arm of the LRPl model is (i) located more 
to the southeast and (ii) more detached from the disk, 
whereas that of the LRP2 model is (i) located more to 
the northwest and (ii) closer to the galaxy center than 
it is observed. The LRP2 gas distribution shows a local 
maximum ~ 1' southeast of the galaxy center. It is not 
excluded that a significant amount of this gas might be in 
molecular form. 

Both models also reproduce qualitatively the overall 
gas kinematics. The two main characteristics of the ob- 
served Hi velocity field, the earlier flattening of the rota- 
tion curve in the southeast and the regular velocity field 
in the extended low surface-density northeastern region 
indicating rotation, are well reproduced. This behavior is 
more pronounced in the LRP2 model snapshot. For a more 
detailed comparison we show position-velocity diagrams 
along the major axis, minor axis, and parallel to the mi- 
nor axis 20" southeast from the galaxy center (Fig. [^T|) . 
Whereas the observed rising rotation curve to the north- 
west is reproduced, the model rotation curve to the south- 
east decreases sightly, whereas the observed Hi rotation 
curve stays constant. The model nicely reproduces the Hi 
data for the cuts along and parallel to the minor axis. The 
LRPl model distribution of the dispersion velocity rea- 
sonably agree with the Hi observations. That of the LRP2 
model shows an asymmetry along the major axis, with a 
larger velocity dispersion to the southeast due to radial 
velocities induced by ram pressure compression, which is 
not observed. 

6.2. The star formation distribution 

The coUisional component of the model consists of 20 000 
gas clouds of different mass that can have inelastic col- 
lisions. For particles of equal size the collision rate, i.e. 
the number of collisions per unit time, is dN/dt — nav 
where n is the particle volume density, a the cross sec- 



tion, and V the particle velocity. The coUisional timescale 
is icoU = {nav)~^ . The mass involved in collisions per unit 
time and unit area is thus S = T,nav, where E is the sur- 
face density. We now assume that the star formation 
rate is proportional to the cloud collision rate cx E. 
In a turbulent selfgravitating gas disk in hydrostatic equi- 
librium n (X fl^ and E oc fi, where ft is the angular ve- 
locity (VoUmer & Beckert 2002,2003). Thus, we obtain 
E* oc E^-^ which is close to the observed Schmidt law (see, 
e.g. Wong & Blitz 2002). To compare the model star for- 
mation rates with Ha observations we gather all collisions 
up to 10 Myr before the timestep of interest. The result- 
ing model star formation rates are shown in Fig. 1221 These 
snapshots can be directly compared to Fig.[71 Both model 
snapshots show a ring of star formation in the inner 5 kpc 
and enhanced star formation along the compressed ridge 
of high column density gas at the southwestern windward 
side. The northwestern Ha spiral arm is less prominent in 
the model. The observed lack of star formation at ~ 1' 
east from the galaxy center is only present in the LRP2 
model. This is directly linked to the discrepancy of the ob- 
served and simulated gas surface density at this location 
discussed in Sec. 16.11 The underlying gas cloud distribu- 
tion of the LRP2 model snapshot is shown in Fig. [531 The 
observed faint outer Ha arm in the west of the galactic 
disk is not present in our model where the star formation 
peaks together with the gas distribution. However, we ob- 
serve a qualitatively similar behavior of the LRP2 model 
star formation distribution southwest from the galaxy cen- 
ter: there are two distinct arms with a clear gap between 
them. The external arm is due to ram pressure compres- 
sion, the internal arm is a spiral arm. Since the Ha emis- 
sion distribution also shows a double arm structure in the 
southwest (Fig. [7]), we suggest that the southern part of 
the compressed Hi ridge consists of two distinct narrow 
gas arms which are smeared out by the 30" resolution of 
the Hi observations. We thus predict that future high res- 
olution Hi observations detect these two arms. 
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Fig. 20. Comparison between the observed Hi and simulated gas moment maps. Upper row: gas distribution; middle 
row: velocity field; lower row: velocity dispersion. 



6.3. The polarized radio-continuum emission 

The comparison between the observed and simulated 
gas surface density distribution, the polarized radio- 
continuum emission, and the sky-projected vectors of the 
large-scale magnetic field is shown in Fig. [241 As explained 
in Sec. 15.31 we lose the large-scale magnetic field in the 
inner disk due to numerical diffusion. Moreover, the ob- 
served polarized emission in the inner disk could be due to 
regular fields or compressed random fields due to dynamo 
action or compressed fields due to density waves, which 
are not part of the model. The model polarized emission is 



due to an enhancement of the large-scale magnetic field by 
shear and radial motions (induction equation) which are 
well present in the model velocity fields. Therefore, while 
the loss of magnetic field energy by diffusion is a numerical 
effect, the gain of magnetic field energy, and thus the en- 
hancement of polarized radio-continuum emission, is real. 
We successfully reproduce the enhanced polarized radio- 
continuum distribution at the edge of the compression re- 
gion. As observed the model polarized emission peaks in 
front of the Hi emission for the LRPl model. However, 
the south-eastern maximum of the LRP2 model is located 
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Fig. 22. Greyscale: model star formation distribution. Contours: model gas distribution. The timestep is indicated in 
the lower right corner of each panel. 



at the backside of the Hi maximum. The magnetic field 
vectors are parallel to this edge for the LRPl, in agree- 
ment with observations, whereas they deviate from this 
direction at the location of the southeastern maximum of 
the LRP2 model. 



7. Discussion 

The first set of simulations with the temporal ram pressure 
profiles used in VoUmer et al. (2001, 2005, 2006) quahta- 
tively reproduce the observed Hi surface density distri- 
bution (LRP of Fig. [m and HRP of Fig. [Ml). However, 
the deduced timescale to peak ram pressure of 120 Myr is 
too short for the galaxy to reach a position near enough 
to the cluster center where the expected ram pressure 
corresponds to the model peak ram pressure. Therefore, 
we made a second set of simulations where we doubled 
the widths of the temporal ram pressure profiles. Since 
the inferred motion of the galaxy within the intraclus- 
ter medium does not permit it to reach a position where 
the peak ram pressure is strong (as for NGC 4438 where 
p^^^ = 5000 cm-3(kms-i)2, HRPl/2), we prefer the 
simulations with moderately strong peak ram pressure 
(p^^^ = 2000 cm-3(kms-i)2, LRPl/2). These dynam- 
ical simulations qualitatively reproduce the observed Hi 
surface density distribution and kinematics (Fig. [20)l. the 
observed massive star formation distribution (traced by 
Ha emission; Fig. [22|l . and the 6 cm polarized radio con- 
tinuum emission distribution (Fig. I24p . 

Within the LRP 1/2 model the current ram pres- 
sure on NGC 4501 is 400 cm-3(kms-i)2 (LRPl model) 
and 800 cm-3(kms-i)2 (LRP2 model). The direction of 
the model galaxy's three dimensional velocity vector is 



(0.48, -0.44, 0.76) (a;-axis: RA, positive to the west; y- 
axis: DEC, positive to the north; z-axis: line-of-sight, 
positive for the far side). Based on this geometry de- 
rived from our simulations, NGC 4501 's orbital velocity is 
1500 kms~^ for both model snapshots LRPl and LRP2. 
A simplified Gunn & Gott criterion (Gunn & Gott 1972) 
yields 

,,2 



■^ISM- 



R 



PiCMWgai 



(4) 



where R is the stripping radius, Sism the Hi column den- 
sity, Vrot the rotation velocity, and Wgas the orbital veloc- 
ity of NGC 4501. This estimate assumes a smooth and 
static intracluster medium. Taking an Hi column den- 
sity of EisM = 10^° cm~^ and a stripping radius of 
12.5 kpc, we derive an intracluster medium density of 
PiCM = 10~^ cm~'^. This is consistent with the intraclus- 
ter medium density at the projected location of NGC 4501 
derived from X-ray observations under the assumption of 
hydrostatic equilibrium picM — 7 10~^ cm~^ (Schindler et 
al. 1999). 

The LRPl model better reproduces (i) the outer edge 
of the Hi distribution compared to the optical isophots 
(Fig. [Tg]) . (ii) the symmetric distribution of the velocity 
dispersion (Fig. [^0)1 . (iii) the distribution of the polarized 
radio-continuum emission (Fig. [24]), and (iv) the galaxy 
has more time to approach the cluster center (320 Myr to 
peak ram pressure). On the other hand, the LRP2 model 
better reproduces (i) the asymmetric Hi distribution along 
the major axis with a larger extent to the southeast, (ii) 
the steepening of the rotation curve to the northwest and 
its flattening to the southeast (Fig. [21]) , and (iii) the gap 
between two arms of star formation at the western edge of 
the galactic disk (Fig. . The timescale to ram pressure 
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Fig. 24. Comparison between the observed and simulated gas surface density distribution (greyscale; resolution: 30"), 
the polarized radio-contimmm emission (contours; resolution: 18"), and the sky-projected vectors of the large-scale 
magnetic field. 



maximum of 190 Myr is smaller than that of the LRPl 
model, but still sufficient to reach a position in the cluster 
where the model peak ram pressure can be reached. Based 
on these arguments we regard the timescale of 320 Myr to 
peak ram pressure as an upper limit. 

Fig. mi shows the linear extrapolation of NGC 4501's 
probable trajectory. With a radial velocity with respect 
to the Virgo cluster mean of 1130 kms^^ a velocity 
within the sky plane of ^ 1000 kms~^ is inferred. Thus, 
the galaxy travels 0.33/0.19 Mpc from its present loca- 
tion within 320/190 Myr. At maximum ram pressure the 
galaxy's trajectory is perpendicular to its distance vector 
to M87. The projected angles between the distance vector 
to M87 and the galaxy's trajectory are 112°/125° and the 
projected distances to M87 are 0.3 Mpc/0.48 Mpc, respec- 
tively. This translates into a three dimensional distance to 
M87 of 0.39 Mpc and 0.54 Mpc at maximum ram pressure, 
respectively. The line-of-sight distance from M87 at this 
moment is then 0.14 Mpc and 0.27 Mpc, respectively. The 
present line-of-sight distance of NGC 4501 is 0.11 Mpc 
in front of M87 for the LRPl model {At = -320 Myr) 
and 0.13 Mpc behind M87 for the LRP2 model {At = 
-190 Myr). 

This overall picture (see also table [3]) makes us con- 
clude that the outer western ridge of star formation is in- 
duced by nearly edge-on ram pressure stripping {i ~ 10°). 
A comparable case is NGC 4654 (Vollmer et al. 2003), 
where nearly edge-on ram pressure stripping (j ~ 0°) leads 
to an enhanced northwestern Ha emission in the compres- 
sion region. This hypothesis can be verified by the deter- 
mination of the star formation history from multiband 
photometry including UV emission and optical spectra. If 
our scenario is correct we expect a recent small starburst 
that created the western outer Ha arm. 

8. Conclusions 

VIVA Hi observations of the Virgo spiral galaxy NGC 4501 
are presented. The outer edge of the Hi distribution fol- 



lows one optical isophot contour except in the northeast, 
where a region of low surface-density gas is discovered. 
Whereas the B band image is asymmetric, the lowest B 
band isophot contour and the H band isophot contours 
are symmetric. The observed southwestern ridge of 6 cm 
polarized radio-continuum emission (Vollmer et al. 2007) 
is located at the outer edge of the Hi distribution. Both 
show the same rapid increase at the edge, which is ex- 
pected in a ram pressure scenario. Along the gas ridge, a 
faint Ha emission ridge is present. The maximum of the 
polarized radio-continuum emission are located within the 
gap between the Ha emission of the spiral arms and the 
faint emission ridge. 

We compare these observational characteristics with 
detailed dynamical models including MHD calculations. 
Models with different galaxy orbits (high and low ram 
pressure) and disk-wind angles between the disk and the 
orbital plane are presented. Model snapshots at differ- 
ent times from peak ram pressure are compared to the 
multi- frequency observations. For the first time we present 
model maps of massive star formation. We qualitatively 
reproduce the gas distribution, kinematics, star formation 
distribution, and polarized radio-continuum emission dis- 
tribution. We conclude that ram pressure stripping can ac- 
count for the main characteristic. The western ridge of en- 
hanced gas surface density and enhanced polarized radio- 
continuum emission is due to ram pressure compression. 
Assuming radial orbits in a static intracluster medium of 
the Virgo cluster (Vollmer et al. 2001), NGC 4501 will 
have its nearest approach to the cluster core, i.e. peak 
ram pressure, in ~ 200 — 300 Myr (Fig. [9]). It is argued 
that the faint western Ha emission ridge is induced by 
nearly edge-on ram pressure stripping. 
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Fig. 21. Comparison between observed Hi (greyscale) and 
simulated position- velocity diagrams for the LRPl (dotted 
line) and LRP2 model (solid line). Upper panel: cut along 
the major axis, northwest is to the right; middle panel: 
cut along the minor axis, northeast is to the left; lower 
panel: cut parallel to the minor axis 20" southeast from 
the galaxy center, northeast is to the left. 



2693/H03/2006/31. We would like to thank the referee for his 
stimulating comments. 



References 

Beck, R., Fletcher, A., Shukurov, A., et al. 2005, A&A, 444, 
739 

Beck, R. 2005, in Cosmic Magnetic Fields, ed. R. Wielebinski 
& R. Beck (Springer: Berlin), p. 41. 



20 



10 



o 

Q 



-1{ 



-20 




-20 -10 

RA (kpc) 

Fig. 23. Model cloud distribution at t 
LRP2 snapshot. 



20 



-190 Myr of the 




Fig. 25. Sketch of the linear extrapolation of NGC 4501 's 
probable trajectory within the Virgo cluster. Contours 
show the Hi surface density of the large spiral galaxies 
around the cluster center (M87; from Cayatte et al. 1990). 
The projected distance between NGC 4501 and M87 is 
2° (0.6 Mpc). The two lines between M87 and the tra- 
jectory correspond to the extrapolated galaxy's projected 
distance at peak ram pressure for the LRP2/HRP2 snap- 
shots (the galaxy is observed 190 Myr before peak ram 
pressure) and the LRPl/HRPl snapshots (the galaxy is 
observed 320 Myr before peak ram pressure). 

Bohringcr, H., Neumann, D.M., Schindler, S., Huchra, J. P. 

1997, ApJ, 485, 439 
Boulares, A. & Cox, D.P. 1990, ApJ, 365, 544 
Cayatte, V., van Gorkom, J.H., Balkowski, C, Kotanyi, C. 

1990, AJ, 100, 604 

Cayatte, V., Kotanyi, C, Balkowski, C, van Gorkom, J.H. 

1994, AJ, 107, 1003 
Chamaraux, P., Balkowski, C, Gerard, E., 1980, A&A, 83, 38 
Chung, A., van Gorkom, J.H., Kenney, J.D.P., & VoUmer, B. 

2007, ApJ, 659, L115 
de Vaucoulcurs G., do Vaucoulcurs A., Corwin H.G., ot al. 

1991, Third Reference Catalogue of Bright Galaxies, (New 



VoUmer et al.: Pre-peak ram pressure stripping in NGC 4501 



21 



York:Springer)(RC3) 
Elmegrccn, B.G. & Falgaron, E. 1996, ApJ, 471, 816 
Elstner, D., Otmianowska-Mazur, K., von Linden, S., Urbanik, 

M. 2000, A&A, 357, 129 
Gavazzi, G., Boselli, A., Donati, A. et al. 2003, A&A, 400, 451 
Giovanelli, R. & Haynes, M.P. 1983, AJ, 88, 881 
Guhathakurta, P., van Gorkom, J.H., Kotanyi, C.G., 

Balkowski, C. 1988, AJ, 96, 851 
Gunn, J.E., Gott, J.R. 1972, ApJ, 176, 1 
Harnett, J., Ehle, M., Fletcher, A., et al. 2004, A&A, 421, 571 
Hummel, E. & Beck, R. 1995, A&A, 303, 691 
Isaaks, E.H., Srivastava, R. 1989, An Introduction to Applied 

Geostatistics, Oxford University Press, Chapter 12. 
Nulsen, P.E.J. 1982, MNRAS, 198, 1007 
Otmianowska-Mazur, K. & VoUmer, B. 2003, A&A, 402, 879 
Schindler, S., BinggeU, B., & Bohringer, H. 1999, A&A, 343, 

420 

Soida, M., Urbanik, M., & Beck, R. 1996, A&A, 312, 409 
Soida, M., Urbanik, M., Beck, R., Wielebinski, R., & 

Balkowski, C. 2001, A&A, 378, 40 
Soida, M., Otmianowska-Mazur, K., Chyzy, K., & VoUmer, B. 

2006, A&A, 458, 727 
Springel, V., Yoshida, N., & White, S.D.M. 2001, NewA, 6, 79 
Stone J.M., Norman M.L., 1992a, ApJS 80, 791 
Stone, J. M., Norman, M. L., 1992b, ApJS, 80, 791 
VoUmer, B., Cayatte, V., Balkowski, C., & Duschl, W.J. 2001, 

ApJ, 561, 708 
VoUmer, B. & Beckert, T. 2002, A&A, 382, 872 
VoUmer, B. & Beckert, T. 2003, A&A, 404, 21 
VoUmer, B. 2003, A&A, 398, 525 

VoUmer, B., Braine, J, Combes, F., Sofue, Y. 2005, A&A, 441, 
473 

VoUmer, B., Soida, M., Otmianowska-Mazur, K., et al. 2006, 

A&A, 453, 883 
VoUmer, B. & Huclitmcicr, H. 2007, A&A, 462, 93 
VoUmer, B., Soida, M., Beck, R. et al. 2007, A&A, 464, L37 
Wiegel W. 1994, Diploma Thesis, University of Heidelberg 
Wong T. & Blitz L. 2002, ApJ, 569, 157 



